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FOREWORD 


Recent  ACIC  activities  in  the  field  of  analytic  photogrammetry  have 
developed  a  requirement  among  technicians  for  knowledge  of  the 
mathematical  essentials  of  analytic  aerotriangulation.  This  manual 
is  one  of  a  series  of  manuals  designed  to  provide  the  mathematical 
knowledge  required  to  understand  current  literature  related  to  ACIC 
activities  (other  manuals  already  printed  include  Technical 
Memorandum  TM-17,  “Determinants  and  Matrices"  and  TM-16 
“Functional  Notation,"  prepared  in  the  Research  Division).  The 
material  presented  reoresents  information  relative  to  the  geometric 
basis  for  solution  to  the  aerotriangulation  problem,  which  is 
fundamentally  the  projective  relations  between  a  point,  a  line,  and 
a  plane  in  space. 

There  are  six  chapters  in  this  manual.  The  first  chapter  gives  a 
geometric  definition  of  the  general  aerotriangulation  problem  and 
reviews  some  fundamental  propositions  from  solid  geometry  which 
are  necessary  for  further  study.  Subsequent  chapters  deal  with  the 
point  in  space,  the  line  in  space,  the  plane  in  space,  and  systems 
of  coordinates  with  matrix  notation  for  transforming  coordinates 
from  one  reference  system  to  another  system.  The  transformation 
from  a  space  coordinate  system  to  a  geodetic  coordinate  system 
based  on  a  spheroidal  earth  is  a  specific  consideration.  The  last 
chapter  deals  with  practical  applications. 


GEOMETRIC  BASIS  OF  ANALYTIC  AEROTRIANGULATION 


1.  INTRODUCTION 

1.1  The  General  Problem. 

The  general  purpose  of  analytic  aerotriangulation  is  to  establish  ground  positions 
from  the  rectangular  coordinates  of  objects  recorded  on  photographs  or  photographic 
plates.  The  aerotriangulation  problem  is  basically  the  same  as  any  other  measuring  prob¬ 
lem  —  precise  measurement  and  the  exact  mathematical  expression  of  the  geometric  rela¬ 
tionships  in  the  problem. 

The  following  sections  are  introductory  to  any  treatment  of  the  analytic  aerotri¬ 
angulation  problem.  The  entire  analytical  solution  is  based  on  a  few  theorems  and  con¬ 
cepts  from  plane  and  solid  analytical  geometry.  This  manual  is  designed  to  assist  the 
reader  in  understanding  the  use  of  these  concepts  in  the  aerotriangulation  problem. 

The  presentation  emphasizes  the  projective  relationships  between  a  point  and 
plane  in  space,  a  line  and  plane  in  space,  and  planes  in  space.  Even  though  these  topics 
are  interrelated  and  could  very  well  be  included  in  one  chapter,  each  topic  (the  point,  the 
line  and  the  plane  in  space)  is  treated  separately  to  show  its  individual  application  to  the 
problem.  Only  those  derivations  necessary  to  visualize  the  geometric  concepts  are  dis¬ 
cussed,  including  the  transformation  of  coordinates  in  space  and  systems  of  coordinates. 

1.2  Fundamental  Definitions  and  Propositions,  The  following  basic  definitions  and 
propositions  from  solid  geometry  are  essential  to  an  understanding  of  the  material  in  the 
in  the  next  four  chapters: 

a.  The  locus  of  a  point  is  the  path  of  the  point  when  it  moves  according  to  some 
definite  law.  In  plane  geometry,  a  locus  is  usually  a  straight  line  or  a  curve.  In  solid 
geometry,  a  locus  is  usually  a  plane  or  a  curved  surface. 

b.  A  straight  line  is  the  locus  of  a  point  in  a  plane  always  moving  in  a  fixed 
direction.  Generally,  a  surface  is  the  locus  of  a  point  in  space  which  satisfies  one  given 
condition.  The  locus  of  a  point  in  space  which  satisfies  two  conditions  is  generally  a 
curve. 

c.  A_ray  is  a  line  extending  an  infinite  distance  in  one  direction  from  a  fixed 

point. 

d.  A  segment  is  a  limited  portion  of  a  line. 

e.  A  plane  is  determined  by: 

(1)  A  straight  line  and  a  point  outside  the  line. 

(2)  Three  points  which  are  not  on  the  same  straight  line. 

(3)  Two  intersecting  lines  or  two  parallel  lines. 

f.  The  orthogonal  projection  of  a  point  on  a  plane  is  the  foot  of  a  perpendicular 
drawn  from  the  point  to  the  plane. 
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g.  The  projection  of  a  line  on  a  plane  is  the  line  containing  the  projections  of  all 
points  of  the  given  line  on  the  given  plane. 

h.  Two  planes  are  parallel  if  they  never  meet,  however  far  they  are  extended. 

i.  One  and  only  one  plane  parallel  to  the  given  plane  can  be  drawn  through  a 
point  outside  a  plane. 

j.  All  the  perpendiculars  to  a  given  line  at  a  point  in  the  line  lie  in  a  plane  perpen¬ 
dicular  to  the  line  at  that  point. 

k.  Two  planes  are  perpendicular  if  they  form  a  right  dihedral  angle. 

l.  If  a  line  is  perpendicular  to  a  plane,  every  plane  containing  the  line  is  perpen¬ 
dicular  to  the  given  plane. 

m.  Through  a  given  line  not  perpendicular  to  a  given  plane,  one  plane  and  only 
one  plane  can  be  drawn  perpendicular  to  the  given  plane. 
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2.  THE  POINT  IN  SPACE 


2.1  Rectangular  Coordinates  in  Space.  The  idea  of  rectangular  coordinates  as  devel¬ 
oped  in  plane  analytical  geometry  may  be  extended  to  space  by  assuming  three  mutually 
perpendicular  planes  intersecting  in  the  lines  XX',  YY',  and  ZZ',  Figure  1. 


The  planes  XY,  YZ  and  XZ  are  coordinate  planes  and  their  lines  of  intersection  establish 
the  y-axis,  z-axis  and  x-axis.  The  positive  directions  are  indicated  by  the  arrowheads  and 
the  point  of  intersection  is  the  origin.  The  eight  sections  of  space  separated  by  the  coor¬ 
dinate  planes  are  called  octants.  Let  P  in  Figure  2  be  in  the  first  octant:  the  octant 
behind  the  first  is  numbered  the  second,  the  octant  to  the  left  of  the  second  is  the  third, 
and  the  octant  in  front  of  the  third  (and  to  the  left  of  the  first)  is  the  fourth.  By  starting 
with  the  octant  directly  below  the  first  and  proceeding  in  a  similar  manner,  the  fifth, 
sixth,  seventh  and  eighth  octants,  respectively,  are  obtained. 
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The  octant  in  which  a  point  is  located  can  be  determined  by  the  signs  or  direction 
of  its  coordinate  axes.  The  following  table  illustrates  the  classification  of  octants  accord¬ 
ing  to  the  signs  of  the  coordinate  axes: 


Octant 

x-axis 

y-axis 

z-axis 

1st 

+ 

+ 

+ 

2nd 

- 

+ 

+ 

3rd 

- 

- 

+ 

4th 

+ 

- 

+ 

5th 

■f 

- 

6th 

- 

+ 

- 

7th 

- 

- 

- 

8th 

+ 

% 


Figure  2 

In  Figure  2,  three  planes  are  drawn  through  the  point  P  parallel  to  the  three  coor¬ 
dinate  planes  and  cutting  the  axes  at  A,  B  and  C.  The  distances,  OA  =  x,  OB  =  y,  and 
OC  “  z  are  the  rectangular  coordinates  of  P.  They  are  equal  numerically  to  the  perpendic¬ 
ular  distances  from  the  coordinate  planes  and  designated  (x,  y,  z)  or  for  a  point  in  the 
fourth  octant  (2,  -3,  4). 
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2.2  Equations  of  a  Point  in  Space,  Since  three  planes  through  the  point  P  in  Figure  2 
are  drawn  parallel  to  the  coordinate  planes,  which  are  mutually  perpendicular,  then 

OA  =  LP,  OB  =  MP  and  OC  -  NP.  The  lines  LP,  MP  and  NP  are  perpendicular  to  the 
ZY,  XZ  and  XY  planes  respectively.  By  letting  LP  =  a,  MP  =  b  and  NP  --  c,  the  three 
equations 

x  =  a,  y  =  b,  z  =  c  (2-1) 

are  the  equations  of  the  point  P.  For  the  point  N,  in  the  XY  plane,  z  =  0  and  the  equa¬ 
tions  of  N  are 

x  =  a,  y=b,  z  =  0  (2-2) 

For  the  point  A,  on  the  x-axis,  both  y  and  z  are  zero,  and  the  equations  of  A  are 

x=a,  y  =  0,z-0  (2-3) 

In  every  case,  three  equations  are  necessary  to  determine  a  point  in  three-dimensional 
space, 

2.3  Distance  Between  Two  Points  in  Space. 
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Let  Pi  and  P2  (Figure  3)  be  two  points  in  space  with  coordinates  Xj,  yj,  Zj  and 
X2,  y2,  Z2  respectively.  By  applying  the  Pythagorean  theorem 

(PiD)^  =  (Pi  A)^  +  (AD)^  =  (x2  -  xj)^  +  (y2  -  yi)^  (2-4) 

and 

(PiP2)'=  (PiD)S  (DP2)'  (2-5) 

where  (DP2)^=  {z^  By  substituting  the  values  of(PjD)^and  (DP^  )^  inequation 

(2-5)  and  denoting  (P^  P^)^  by  d^,  the  distance  between  two  points  in  space  is  determined 
by  the  formula 

d^  ==  (X2  -  Xi)2  +  (y2  -  yi)2  +  (Z2  -  Zi)2 

d  ===  V(X2  -  Xj)^  +  (y2  -  YiV  +  (Z2  -  Zi)2  (2-6) 

If  Pi  is  the  origin  (0,0,0)  and  the  coordinates  of  are  x,y,z,  then 

d  =  V(x  -  0)2  +  (y  -  0)"  +  (z  -  0)2  or 

d  =  Vx2  +  y2  +  z2  (2-7) 

2.4  Orthogonal  Projection  of  a  Point  in  Space.  The  orthogonal  projection  of  a  point 
on  a  coordinate  plane  is  the  foot  of  the  perpendicular  from  the  point  to  the  plane. 


Figure  4 

In  Figure  4,  L  is  the  orthogonal  projection  of  the  point  P  upon  the  XZ  plane, 

M  is  the  orthogonal  projection  of  P  upon  the  XY  plane,  and  N  is  the  orthogonal  projection 
of  P  upon  the  YZ  plane. 
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3.  THE  LINE  IN  SPACE. 


3.1  Projections  of  a  Line  Upon  the  Coordinate  Planes.  Two  non-parallel  planes  in 
in  space  will  intersect  in  one  and  only  one  straight  line.  The  projection  of  a  line  Pj  , 
upon  any  one  of  the  three  coordinate  planes  is  the  line  of  intersection  of  the  coordinate 
plane  with  a  plane  containing  the  line  P^  Pg .  The  projection  is  the  orthogonal  projection 
of  the  line  P|  P2  if  the  projecting  plane  is  perpendicular  to  the  coordinate  plane  (Fig¬ 
ure  5). If  the  projecting  plane  is  not  perpendicular,  the  projection  is  oblique. 

% 


3.2  Direction  Numbers  and  Direction  Cosines  of  a  Line.  Even  though  two  lines  in 
space  do  not  intersect,  it  is  convenient  to  define  the  angle  between  them  as  the  angle 
between  intersecting  lines  which  are  respectively  parallel  to  them.  Let  O’P* Figure  6  be 
any  directed  line  in  space,  and  OP  be  a  line  through  the  origin  0  which  has  the  same 
direction  and  is  parallel  to  O'  P'.  The  angles  a  ,  and  y  are  arigles  which  OP  makes  with 
the  coordinate  axes.  These  are  direction  angles  and  their  cosines  are  direction  cosines 
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Figure  6 

of  the  line  OP.  The  latter  will  be  denoted  by  \ ,  |ut ,  and  v  respectively.  Let  P  be  denoted 
by  the  coordinates  (a,  b,  c),  then  in  the  triangle  OPA 

A  =  cos  a  =  (3-1) 

Similarly  in  triangles  OBP  and  OCP  respectively, 

u  =  cos  6  =  V  !=  cos  y  =  .  (3-2) 

OP  OP 

Since  OP  is  the  diagonal  of  a  parallelepiped  whose  edges  are  a,  b  and  c, 

OP  =  Va2  +  +  c2.  (3-3) 
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By  substituting  the  value  of  OP  in  Equations  (3-1)  and  (3-2),  the  direction  cosines  can 
be  expressed  as 

A  «  cos  a  _ 2 _ 

+  Va^  +  b^  +  c^ 

fi  =  cos  P  =  _ h _  (3-4) 

+  Va^  +  b2  ^  q2 

V  =  cos  y  =  _ c _ 

±  Va2  +  b2  +  c2 

By  squaring  each  member  of  the  latter  equations  and  adding  the  results,  the  fundamental 
formula  of  space  trigonometry  is  obtained: 

cos  ^  a  -f  cos"  p  +  cos  y  =  1.  or  A^  +  +  1^2  =  1  (3-5) 

The  sum  of  the  squares  of  the  direction  cosines  of  a  line  is  equal  to  unity. 

Any  three  numbers  proportional  to  the  direction  cosines  of  a  line  are  the  direction 
numbers  of  the  line.  By  definition 

cos  a  _  cos  p  _  cos  y  _  j. 
a  b  c 

If  r  represents  the  common  value  of  the  ratios,  then, 

cos  a  =  ar,  cos  p  =  br,  cos  y  =  cr.  (3-6) 

Squaring,  adding,  and  applying  Equation  (3-5) 

(a^  +  b^  +  c^)  =  1,  and 

r  =  1 

+  Va^  +  b^  +  c^ 

Substituting  this  value  for  r  in  Equation  (3-6),  formulas  for  direction  cosines  (3-4)  are 
obtained.  (The  importance  of  this  derivation  is  that  any  three  numbers  a,  b,  and  c  deter¬ 
mine  the  direction  of  a  line  in  space.  This  direction  is  the  same  as  that  of  the  line  joining 
the  origin  and  the  point  (a,  b,  c)  when  the  positive  sign  of  the  radical  is  chosen,  other¬ 
wise  the  direction  is  reversed.) 

If  one  direction  number  is  zero,  the  line  is  perpendicular  to  the  corresponding  axis 
of  the  coordinates.  If  two  are  zero,  it  is  parallel  to  the  remaining  axis.  The  following 
considerations  will  guide  the  determination  of  signs  for  the  radical  in  Equation  (3-4): 

a.  If  a  line  cuts  the  XY-plane,  it  will  be  directed  upward  or  downward  according 
as  cos  y  is  positive  or  negative. 
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b.  If  a  line  is  parallel  to  the  XY-plane,  cos  y  =  0,  and  the  line  will  be  directed 
in  front  or  in  the  back  of  the  ZX-plane  according  as  cos  a  is  positive  or  negative. 

c.  If  a  line  is  parallel  to  the  x-axis,  cos  =  cos  y  =  0,  and  its  positive  direc¬ 

tion  will  agree  or  disagree  with  that  of  the  x-axis  according  as  cos  a  =  1  or  -1. 

3.3  Equation  of  a  Line  in  Space.  In  plane  analytic  geometry,  the  equation  of  a 
straight  line  on  the  XY  plane  is: 

Ax  +  By  +  C  =  0,  (3-7) 

where  A,  B,  and  C  are  constants. 


Figure  7 

C 

Let  the  line  PQ  intersect  the  y-axis  at  the  point  y  ~ - ,  and  makes  an  angle  6  with 

B 

the  x-axis  such  that 


tan 


B  ^ 


Figure  7. 


(3-8) 


Consider  this  same  equation  in  three-dimensional  space.  The  same  relationship 
between  x  and  y  holds  with  no  limit  on  the  z-coordinate.  The  locus  of  the  equation  is  a 
plane  perpendicular  to  the  XY  coordinate  plane  and  makes  the  angle  9  with  the  XZ  coor¬ 
dinate  plane  (Equation  3-8,  and  Figure  8). 
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Figure  8 

The  theorem  from  solid  geometry  which  states  that  two  non-parallel  planes  inter¬ 
sect  in  a  straight  line  leads  to  the  conclusion  that  one  form  of  the  equations  of  a  straight 
line  can  be  developed  from  any  two  of  the  three  projections  of  the  line  on  the  coordinate 
planes.  Thus  any  two  of  the  following  three  equations  are  the  equations  of  a  line  in  space: 

AjX  +  B^y  +  Cj  =  0 

A2X  +  B2Z  +  C2  =  0  (3-9) 

A^y  +  B^  z  +  Cg  —  0 
where  A,  B,  and  C  are  constants. 
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In  Figure  9,  let  Pj  and  Pg  be  two  known  points  on  a  line  of  known  length  d  with 
known  direction  cosines  a  ,  P ,  and  y  .  Let  Q  be  a  variable  point  with  coordinates  x,  y,  z. 

Then 

^  ~  =  cos  a  , 

d 


y  “Yi 
d 


cos  p  , 


(340) 


^  ~  =  cos  y  , 

d 

and  by  simple  transformation 
X  =  Xi  +  d  cos  a 

y  =  yj  +  d  cos  P  (3-11) 

z  =  Zi  +  d  cos  y 

The  equations  in  (3-11)  are  the  parametric  form  of  the  equations  of  a  line  in  space, 
solving  Equations  (3-11),  for  d  and  equating  the  results: 

^  =  y  -  yi  =  z  (3.12) 

cos  a  cos  p  cos  y 

which  yields  the  symmetric  form  of  the  equations  of  a  line  in  space.  If  a,  b,  and  c  are 
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direction  numbers  the  symmetric  form  may  be  written  in  the  form 


X  ^  y-yi  ^  z  -z^  , 

a  b  c 


(3-13) 


since  direction  numbers  are  proportional  to  the  direction  cosines. 

The  equation  of  the  straight  line  passing  through  two  known  points  Pi  and  P2  is 
X  -Xi  ^  y  -yi  =  z  -Zi  (3.14) 

X2  “Xi  y^  -yi  Z2  -Zj 

Since  the  denominators  in  Equation  (3-14)  are  direction  numbers,  the  two-point  form  fol¬ 
lows  directly  from  the  symetric  form. 

The  simultaneous  equations 

X  =  mz  +  p,  y  =  nz  +  q  (3-15) 

define  a  line  as  the  intersections  of  its  projecting  planes  on  XZ  and  YZ,  This  is  the 
projection  form  of  the  equations  of  a  line.  Separately,  each  equation  is  the  equation  of  a 
projecting  plane.  The  equations  are  identical  to  the  point  slope  form  of  the  equation  of  a 
line  on  a  plane,  where  m  and  n  are  slopes;  i.e.,  tan  ^1  =  m,  tan  62  =  n,  and  p  and  q  are 
intercepts.  Solving  for  z  and  equating  results 

X  ~P  =  y  -q  =  (3-16) 

m  n  1 

Comparison  with  the  symmetric  form  shows  that  in  the  projection  form  the  line  passes 
through  (p,  q,  o)  and  m,  n  and  1  are  direction  numbers.  Since  m,  n,  and  1  are  direction 
numbers  and  proportional  to  the  direction  cosines,  the  relationship  of  the  direction  cosines 
to  the  slopes  of  the  projecting  planes  can  be  derived  directly  from  the  equations.  Let  k 
be  the  constant  of  proportionality  and, 

mk  =  cos  a,  nk  =  cos  yff,  k  =  cos  y  (3-17) 

Squaring  and  adding 


m^k^  +  n^k^  +  k^  =  k2(m2  +  n^  +1) 

=  cos^  a  +  cos2  +  cos2y  =  1 


(3-18) 


Solving  for  k 


m  -  +  n^  +  1 


Vm2  +  n^  +  T 


(3-19) 
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Substituting  this  value  of  k  in  Equation  (3-17) 


m 

cos  a  =  ■■ 

Vm2  +  n2  +  1 


n  ** 

COS  P  ~  — 

Vm2  +  n2  +  1 

1 

cos  y  =  -— = 
V  +  n2  +  1 


(.3-20) 


3.4 


and 


Condition  for  Intersection  of  Two  Lines  in  Space. 

Let 

X  =  mz  +  p 

(.3-21) 

y  =  nz  +  q 

(.3-22) 

X  =  m'z  +  p' 

(3-23) 

y  =  n'z  +  q' 

(.3-24) 

be  the  equations  of  any  two  lines  in  space.  In  general,  these  lines  will  not  intersect, 
since  the  coefficient  m,  n,  p  and  q  are  entirely  arbitrary.  The  four  equations  contain  three 
unknowns,  x,  y,  z  which  are  not  always  consistent,  even  though  the  lines  in  space  do  not 
intersect,  the  respective  projection  will,  unless  they  are  parallel.  By  eliminating  x  in 
Equations  (3-21)  and  (3-23),  and  y  in  Equations  (3-22)  and  (.3-24),  then, 


and 


(3-26) 


Equations  (3-25)  and  (3-26)  give  the  value  of  the  z-coordinate  of  the  point  of  intersection 
of  the  xz-  and  yz-  projections  of  the  two  lines.  These  values  usually  are  unequal;  how¬ 
ever  if  the  lines  intersect,  the  values  of  z  will  be  equal.  Conversely,  if  the  values  are 
equal,  the  lines  intersect.  The  necessary*  and  sufficient  condition  that  two  lines  in  space 
shall  intersect  is  determined  by  equating  (.3-25)  and  (3-26),  such  that 

(3-27) 

m  -  m'  n  n' 

Equation  (3-27)  gives  the  condition  for  intersection  and,  at  the  same  time,  the 
value  of  the  z-coordinate  of  the  point  of  intersection.  To  find  the  other  coordinates, 
substitute  the  value  of  z  from  (3-25)  in  (3-21)  or  (3-23)  and  from  (3-26)  in  (3-22)  or  (3-24), 
the  resulting  values  will  be, 
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X 


mp*-  m'p 
m  -  m* 


(3-28) 


nq*  -  n»q 
n  -  n» 


(3-29) 


then, 


As  a  numerical  example  of  the  material  presented  in  this  section,  let 
x=3z-2  x  =  -7z+8 

c 

y  =  -4z  +5  y  =  z 


m  =  3,  p  =  -2,  n  =  -4,  q  =  5 

= -7,  p'  8,  n*  =  1,  q»  =  0 

Substituting  these  values  in  (3-27), 

8  4-2  0^5 

3+7  “-4-1 


or  1=1 

The  example  shows  that  the  lines  do  intersect  at  the  point  where  the  z-coordinate  is  L 
By  substituting  the  value  of  1  in  either  set  of  xz-  and  yz  projection  equations,  x  =  1  and 
y  =  1  and  the  two  lines  intersect  at  the  point  (1,  1,  1). 


3.5  The  Angle  Between  Two  Lines. 
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Let  01  ^  and  OP 2  (-see  Figure  10)  be  two  directed  lines  in  space,  such  that 

(OPj)'  x,2  +yi2  +Zi2 

and  (OP2)^  =  X2^  ^  y2^  +  Z2^  i 

also  (PiP2^^  =  (Xj  *-  X2)^  +  (Yi  -  y2)^  -Z2)^. 

Applying  the  Law  of  Cosines  from  Plane  Trigonometry, 

(PiPjj^  =  (OPi)^  WOPj)'  -2  (OPp  (OP2)  cos  P,  OP2 

(OPi)'  +  (OP2)' 

where  cos  Pj  0  P2  =  - • 

2  (OPii  (OP2) 


(3-30; 


(;3-31) 

(3-32) 


By  substituting  in  (3-31)  the  values  obtained  in  (3-30)  and  letting  angle  P,  OP2  be 
represented  by  6,  then; 


cos  e  =  ■  (.3-33) 

(OP,)  (OP2) 

From  the  definition  discussed  with  direction  cosines: 

Xi  =  OPi  cos  ai,  Yi  =  OPi  cos  Pi  ,  =  OPi  cosyi 

and  X2  =  OP2  cos  q2  ,  >'2  =  OP2  cos  ^2 »  ^2  =  0^2  cos  72- 

Substituting  these  values  in  (3-33),  the  following  formula  is  obtained  for  the  angle  between 
two  lines  in  terms  of  the  direction  cosines  of  the  lines: 

cos  0  =  cosaj  cos  a2  cos cos  ^2*  cos  yj  cos  72 

or  cos  ^  M2  ^1^2 

For  parallel  lines: 

cos  =  cos  02 »  cos  =  cos^2»  cos  yj  =  cos  72 
or  Xj  =  X2,  gi  =  g2»  =  ^'2 

For  perpendicular  lines,  Q  =  90°  and  cos  ^  =  0  : 

COSO]  cos  02  “  cos  Pi  cos  /32  -^COSyi  COS  72  =  0* 

or  Xi\2  -p  gig2  ^1^2  =  O' 
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4.  THE  PLANE  IN  SPACE 


4.1  Definitions.  A  plane  is  defined  as  a  surface  such  that  a  line  connecting  any  two 
points  of  it  will  lie  wholly  within  the  surface.  A  plane  may  be  defined  as  the  surface 
generated  by  a  line  rotating  about  another  line  to  which  it  is  perpendicular.  The  rotating 
line  is  called  the  generator  and  the  fixed  line,  or  axis,  the  director.  A  plane  may  be  deter¬ 
mined  by  three  points  (not  on  the  same  line),  or  by  a  line  and  an  external  point,  or  by  two 
intersecting  or  parallel  lines. 

4.2  Equations  of  a  Plane, 

4,2.1  Normal  Form.  The  normal  form  of  the  equation  of  the  plane  is 
X  cos  a  +  y  cos  p  +  z  cos  y  =  p 

This  important  form  gives  the  relationship  between  the  coordinates  of  all  points  in  a 
plane  in  terms  of  the  perpendicular  distance  from  the  origin  to  the  plane,  and  the  direction 
cosines  of  this  perpendicular.  Figures  11a  and  11b. 

In  the  diagram  (11a)  the  plane  LMN  represents  the  position  of  the  plane  in 
one  of  the  octants.  The  line  OD  is  perpendicular  to  the  plane  and  intersects  the  plane 
at  D.  The  length  of  the  normal  OD  is  denoted  by  p  which  is  considered  a  positive  quan¬ 
tity.  The  angles  a ,  p  and  y  are  the  direction  angles  of  the  plane  director,  p  .  When  the 
value  of  p  and  the  direction  cosines  are  known,  the  point  D  can  be  located,  and  the  plane 
is 


% 
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Figure  lib 

constructed  as  a  plane  through  D  and  perpendicular  to  OD.  P  is  any  point  in  the  plane 
LMN,  whose  coordinates  are  x,  y  and  z.  The  triangle  POD  is  always  a  right  triangle 
regardless  of  what  point  in  the  plane  is  chosen  as  P. 

,  If  one  of  the  direction  angles  is  a  right  angle,  p  is  perpendicular  to  one  of 
the  coordinate  axes  and  the  plane  is  parallel  to  that  axis.  For  example,  if  y  =  90®  the 
equation  becomes 

X  cos  a  +  y  cos  )3  =  p 

(the  equation  of  a  plane  parallel  to  the  Z-axisj.  If  two  of  the  direction  angles  are  right 
angles  the  equation  reduces  to  the  form 

x  cos  a  =  p  ,  y  cos  p  =  p  ,  or  z  cos  y  =  p  . 

If  p  and  y  are  right  angles,  p  coincides  with  the  x-axis  and  the  plane  is  perpendicular  to 
the  x-axis. 

4.2.2  General  Form.  Every  equation  of  the  first  degree  can  be  written  in  the 

form: 

Ax  +  By  Cz  -4-  D  =  0 

and  any  equation  of  this  form  is  which  A,  B,  C  are  not  all  zero  represents  a  plane.  This 
general  equation  of  a  plane  is  reduced  to  the  normal  form  by  the  following  rules: 

Rule  I.  Change  all  signs  so  that  the  absolute  term  of  the  given  equation  is 
negative  in  the  first  member,  (p  is  always  considered  positive). 

Rule  2.  Divide  by  the  square  root  of  the  sum  of  the  squares  of  the  coeff'icients  of 
coefficients  of  x,  y,  z. 
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The  normal  form  of  the  general  equation  is 


( 


Ax  +  By  +  Cj  +  D 


) 


=  0 


Va^  +  b"  + 

The  value  of  the  direction  cosines  of  the  plane  director  are; 


cos  a  -  + 


wTb^Tc^ 


cos  fi 


B 


and 


r2 


cos  y  “  + 


C 

VA2  4-  B2  +  C2 


•  The  value  of 


P 


^  D 
=  "  VFTb^T^ 


The  upper  sign  is  to  be  used  when  the  absolute  term  of  the  given  general  equation  is 
positive  and  the  lower  sign  when  it  is  negative. 


Several  special  cases  of  the  general  equation  are  listed  below: 


Equation 

A,  +  By  +  D  =  0 

A,.C,^D  =0 
By  -f  C,  -  D  =  0 
A^,  -f  D  =  0 

By  -h  D  0 

C,  +  D  =  0 


Conditions 


±  xy  plane 

ii 

z-axis 

i_  XZ  plane 

II 

y-axis 

J_  yz  plane 

Ii 

x-axis 

J_  x-axis 

Ii 

yz  plane 

±  y-axis 

1! 

XZ  plane 

J_  z-axis 

Ii 

xy  plane 

4.2.3  Intercepts  and  Traces.  The  intercepts  of  the  plane  on  the  respective  axes 
are  the  distances  from  the  origin,  0,  to  the  point  of  intersection  L,  M,  and  N  of  the  plane 
with  the  coordinate  axes  x,  y,  z.  The  traces  of  the  given  plane  are  the  lines  in  which  a 
plane  intersects  the  coordinate  planes.  The  lines  LM,  MN,  and  NL  (Figure  llaj  are  the 
traces  of  the  ;.'ane  on  the  coordinate  planes  XY,  YZ,  and  XZ. 

The  symmetric  (intercept)  form  of  the  general  equation  is 


L..  *  y  z  ^  Y 
D  _  D  p 
A  B  C 


•  u  D  D  j  D  .  , 

1  he  intercepts  on  the  x,  v,  z  axes  are  -  T’»  “  5"  “  r  respectiveiv. 

A  D  L 

The  slope  form  of  the  general  equation  is: 


A 


B 


z  .  =  0 
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Hfince,  the  general  equation  of  the  first  degree  represents  a  plane  where  the  slopes  of  the 
xz  and  yz  projections  of  its  plane  director  are  A/C  and  B/C  respectively  and  the  intercept 

on  the  z-axis  is  equal  to  -  For  example,  in  the  plane  2x  -  4y  -  3z  4-  12  =  0  to  obtain 
the  symmetric  form,  transpose  the  absolute  term  to  the  second  member  and  divide  by  the 
transposed  absolute  term  to  obtain  ^  ^  Hence,  the  intercepts  are  x  =  -6, 

y  =  3,  z  =  4.  The  plane  is  shown  in  Figure  12. 

In  the  same  example,  to  obtain  the  slope  form,  divide  the  given  equation  by 
the  coefficient  of  z  and  obtain: 

-?.x  +  ly  +  Z-  4  =  0 
3  3 

The  slopes  of  the  xy  and  yz  projections  of  its  plane  director  are  -  2/3  and  4/3  respec¬ 
tively  and  the  intercept  on  the  z-axis  is  4. 


% 


The  plane  4x  4-  5y  -  20  =0  (Figure  13)  has  the  intercepts  x  =  5  and  y  =  4  but  no  inter¬ 
cept  on  the  z-axis  because  it  is  parallel  to  that  axis.  The  traces  are  4x  5y  -  20  =  0, 
X  =  5,  y  =  4. 
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Figure  13 

The  symmetric  or  intercept  form  of  the  equation  of  a  plane  is  the  roost  use¬ 
ful  form  to  plot  the  position  of  the  plane.  The  values  of  the  intercept  are  established  and 
plotted  on  the  respective  axes.  Connecting  the  ends  of  the  intercepts  gives  the  traces  of 
the  plane.  In  the  plane  2x  -  4y  -  3z  +  12  =  0  the  trace  in  the  XY  plane  (z  =  0)  is 
2x  -  4y  +  12  =  0.  Similarly,  the  trace  in  the  XZ  plane  is  2x  -  3y  +  12  =  0  and  the  trace 
in  the  YZ  plane  is  4y  +  3z  -  12  =  0.  The  slope  form  of  the  normal  equation  is 


cos  a  cos  B  p 

- X  +  - y  +z  =  - 

cos  7  cos  7  COS  7 

where  the  coefficients  of  x  and  y  are  the  slopes  of  the  plane  director. 


To  change  the  general  Equation  2x  -  4y  -  3z  +  12  =0  to  the  normal  form, 
apply  Rule  1  and  change  all  signs  because  the  sign  of  the  absolute  term  in  the  left  mem¬ 
ber  is  positive  and  obtain  -  2x  *  4y  +  3z  -  12  =  0.  To  apply  Rule  2,  divide  by 


V29'and  get  the  desired  normal  form, 


4 

-  X  +  — — r  y  + 

\f^  \/2r 


z 


so  that 


cos  a 


=  “  — cos  B 


4 


=  — -  cos  y 

\  29 


3  ,  12 

— - and  p  =  - 

\  29 
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4.3  Distance  Between  a  Point  and  a  Plane.  The  distance  from  any  given  plane  to  any 
point  in  space  is  found  by  reducing  the  equation  of  the  given  plane  to  the  normal  form, 
and  substituting  in  the  normal  form  the  coordinates  of  the  given  point  for  x,  y  and  z.  The 
resulting  numerical  quantity  expresses  the  required  distance  (Figure  14).  To  find  the  dis¬ 
tance  of  the  point  (-3,  4,  -5)  from  the  plane  x  +  2y  -  3z  +  8  =  0  reduce  the  equation  of  the 
plane  to  normal  form  (Section  4.2.2J  and  obtain 

~  X  2y  -f  3z  -  8  _  Q 

vTJ 

Substituting  the  coordinates  of  the  given  point  x,  y  and  z  in  the  normal  form  obtain 
+  .3  -  8  -  15  -  8 

Vli 

distance  =  -  or  -  2  Vl4 

vTI 

The  result  means  that  the  actual  distance  from  the  given  plane  to  the  point  is  2  \  14  units 
to  the  scale  of  the  drawing  and  that  the  point  is  on  the  same  side  of  the  given  plane  as 
the  origin.  If  the  given  point  and  the  origin  are  on  opposite  sides  of  the  given  plane,  the 
result  is  positive  because  p  is  considered  a  positive  quantity. 

In  Figure  14,  P  (x',  y',  z’j  is  any  given  point  on  plane  L'M'N'  vrhich  is  parallel 
to  the  given  plane  LMN.  The  line  OD  =  p  is  the  normal  to  the  given  plane  LMN  from  the 
origin  0.  The  resulting  line  OD  intersects  the  plane  L'M’N’  in  D’.  OD  ’  =  p' .  The  line 
PK  is  the  perpendicular  distance  from  P  to  the  given  plane  and  the  required  distance  from 
the  given  plane  to  the  point  P. 


KP  =  DD’  =  OD’-OD  = 


p  -  p 


Figure  14 


Since  the  second  plane  L'M’N'  passes  through  the  given  point  P,  (x',  y'  ,  z*)  its 
coordinates  satisfy  its  equation, 

giving 

X  ’  cos  a  +  y '  ccs  P  +  cos  y  =  p ' . 

Substituting  this  value  of  p*  in  KP  =  p'  -  p 

KP  =  X*  cos  Q  +  y '  cos p  +  z'  cos  y  -  p  . 

This  equation  is  the  same  as  the  normal  form  of  the  equation  of  the  given  plane  with  the 
coordinates  x'  y'  and  z*  substituted  for  x,  y,  and  z. 


4.4  Angles  between  Two  Planes.  The  angle  between  two  planes  is  the  same  as  the 
angle  between  the  directors  of  the  two  planes >  Figure  15.  The  formula  for  finding  the 
angle  0  between  the  two  planes  whose  general  equations  are 

+  By  +  +  D  =:  0  and  A^ x  +  B, y  +  Cj z  +  D  *  0  is 

AA,  BB,  +  CC' 

cos  e  =  — ^ — ~=L= 

Va^  +  Va^'^  4-  B,'  - 

To  find  the  angle  between  the  tv,o  planes  -  x  +  7y  =  11  and  3x  -  4y  -  5z  ==  10 
substitute  in  the  formula 

(^1)  (3)  +  (7)  f4j  J  (OJ  (o)  25 

cos  d  = - —  --  - or  ;^2 

-f  7  ^  +  0  ^  Vs  ^  +  4  5  ^  50 

and  the  angle  is  60"^. 


1 


The  angle  between  two  planes  is  the  same  as  the  angle  between  the  respective  normals, 
ODj  and  OD2.  The  planes  are  not  shown  in  the  drawing. 

4.5  Parallel  and  Perpendicular  Planes.  Two  planes  whose  equations  are 

Ax  +  Bvn-Cz-r-D  =  0 

A’x  +  B’y  -  C’z  -  D’=  0 

are  parallel  when  and  only  when  the  coefficients  of  x,  y  and  z  are  proportional. 

^  .  A  B  C 

that  IS - = - = - 

A’  B’  C’ 
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For  example  the  two  planes 
2x  +  3y  -  z  =0 
and  4x  +  6y  -  2z  +  8  =  0 
are  parallel  because 

2^3  1 

4  6  2 

Two  planes  are  perpendicular  when  and  only  when  the  sum  of  the  products  of  the 
coefficients  of  x,  y,  and  z  is  zero,  that  is 

AA^  f  BB'+CC*  -  0 

For  example  the  two  planes 

2x  +  3y  -  z  =  0 

and  3x-y  +  3z+2  =  0 

are  perpendicular  because 

(2)(3)+(3)(-l)+U)(3)  =  0 

4.6  Special  Planes.  A  plane  whose  equation  has  the  form 

Ax  +  By  +  D  =  0  is  perpendicular  to  the  XY-plane 

By  +  Cz  +  D  =  0  is  perpendicular  to  the  YZ-plane 

Ax  +  Cz  +  D  =  0  is  perpendicular  to  the  ZX-plane. 

That  is,  if  one  variable  is  lacking,  the  plane  is  perpendicular  to  the  coordinate  plane 
corresponding  to  the  two  remaining  variables. 

A  plane  whose  equation  has  the  form 

Ax  +  D  =  0  is  perpendicular  to  the  axis  of  x 

By  +  D  =  0  is  perpendicular  to  the  axis  of  y 

Cz  +  D  =  0  is  perpendicular  to  the  axis  of  z. 

That  is,  with  only  one  variable,  the  plane  is  perpendicular  to  that  axis.  For  example,  the 
plane  Ax  +  D  =  0  is  perpendicular  to  both  the  XY-plane  and  ZX-plane,  and  hence  is  also 
perpendicular  to  the  line  of  their  intersection. 


5.  THE  COORDINATE  TRANSFORMATION  IN  SPACE  AND  SYSTEMS  OF  COORDINATES 


5.1  Transformation  of  Rectangular  Coordinates 

5.1.1  Translation  of  the  Axes.  Given  two  parallel  sets  of  coordinate  axes,  XYZ 
system  and  X*Y*Z* system,  Figure  16,  with  the  point  P  common  to  both  coordinate  sys¬ 
tems.  The  origin  O’  of  the  X' Y'Z*  system  is  the  point  ,  yo,  z^)  of  the  XYZ  system,  and 
the  origin  0  of  the  XYZ  system  is  the  point  (-  x^,,  -  yo,  -  Zq)  of  the  X'Y’Z ’system.  As  in 
plane  analytical  geometry,  the  following  equations  are  obtained: 

x’  X  -  Xq  X  =:  x'  +  Xo 

y'  =  y  -  yo  y  -  y'  +  yo  (5-1) 

Z'=Z“Zo  Z^Z'+Zq 


Figure  16 

If  the  expressions  for  x,  y,  and  z  in  the  second  set  are  substituted  in  an 
equation  in  x,  y  and  z,  the  corresponding  equation  in  x*,  y',  z‘is  obtained.  A  transforma¬ 
tion  such  as  (5-1)  is  the  translation  of  axes. 

5.1.2  Rotation  of  Axes.  Let  the  coordinates  of  a  point  P,  referred  to  a  set  of 
rectangular  axes  OX,  OY,  OZ,  be  x,  y,  z,  and  referred  to  another  system  OX’,  OY',  OZ' 
having  the  same  origin,  be  x\  y',  z*.  Let  x'  =  OL’,  y'  ==  L'M',  z'  =  M’P  (Figure  17);  and 
let  the  direction  cosines  of  Ox' ,  referred  to  Ox,  Oy,  Oz,  be  u  ^  ;  those  of  OY'be 

»  ^2  »  ^3  >  M3  »  **3  »  "here 
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X .  =  cos  a . 


/Xj  ^cosp. 

V.  =  cos  y.  (i  =  I,  2,  3). 
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The  projection  of  OP  on  the  OX  axis  is  x.  The  sum  of  the  projections  of  OL’, 
L'M’and  M'P  is  Xj  x’  +  X2y'  +  X^z'.  This  follows  from  the  theorem  that  the  length  of  the 
projection  of  a  segment  of  a  directed  line  on  a  second  directed  line  is  equal  to  the  length 
of  the  given  segment  multiplied  by  the  cosine  of  the  angle  between  the  lines.  The  second 
and  third  equations  are  obtained  in  a  similar  manner 

X  =  X^x’-H  X2  y’  +X3  z‘ 

y  =  /x^x'  +  /X2y'  +  /X3Z' 

z  =  V  +  V 

where  X.  ,  /x.  ,  v.  are  the  direction  cosines  of  OX,  OY,  and  OZ  with  respect  to  the  OX', 
OY',  and  OZ’  axes  respectively.  Project  OP  and  OL  =  x,  LM  =  y  and  MP  =  z  on  OX',  OY', 
and  OZ then 

X'  =  ^iX  + 

y '  =  X2X  fi^y  *  z  (5-3) 

Z’  =  X3X  -f  .1X3}*  ^3  z 

O- 


The  systems  of  equations  (5-2)  and  (5-3)  are  expressed  conveniently  by  the 
following  diagram 


X* 

y* 

Z» 

X 

X2 

X3 

y 

Ml 

z 

•-I 

^2 

>^3 

5,2  Systems  of  Coordinates. 


5.2.1  Polar  or  Spherical  Coordinate  Systems, 
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Let  P’  be  an  orthogonal  projection  of  P  on  the  XY-plane  as  shown  in 
Figure  20.  The  position  of  P  is  defined  by  the  distance  r,  the  angle  <>  =  ZOP  which  the 
line  OP  makes  with  the  Z-axis,  and  the  angle  B  (measured  by  the  angle  XOP')  which  the 
plane  through  P  and  the  Z-axis  makes  with  the  plane  XOZ.  The  segment  OM  is  the  x 
distance  along  the  OX-axis,  ON  the  y  distance  along  the  OY-axis  and  OS  the  z  distance 
along  the  OZ-axis.  The  numbers  r,  <> ,  0  are  called  the  spherical  coordinates  of  P;  some 
writers  call  them  polar  coordinates  in  space.  The  length  r  is  the  radius  vector,  the  angle 
<l>  is  the  co-latitude,  and  6  is  called  the  longitude. 
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5.2.2  Transformation  from  Spherical  Coordinates  to  Rectangular  Coordinates  and 
Vice  Versa. 


If  P  =  (x,  y,  z),  then  from  Figure  18, 

OP’  =  r  cos  (90  -  =  r  sin 

In  triangle  OMP*,  x  =  OM  =  OP’  cos  B  =  r  sin  cos  B  . 

In  triangle  ONP*,  y  =  ON  =  OP’  sin  B  r  sin  sin  B . 

In  triangle  OP’P,  z  =  PP'  =  OP '  cos  =  r  cos  ,  or  for  rectantular  coor¬ 
dinates 

X  =  r  sin  <f>  cos  B 
y  =  r  sin  <f>  sin  B 
z  =  r  cos  <l> 

On  solving  for  the  spherical  coordinates  ,  B ,  r  in  terms  of  the  rectangular 
coordinates  x,  y,  z 

r  =  Vx^  +  y  2  +  z2  (Distance  between  two  points 
in  space) 

<p  =  arc  cos  _ 5 _ 

Vx2  y2  +  z2 

B  =  arc  tan  x/y 

5.2.3  Definition  ■  Right-Handed-Rectangular  Coordinate  System. 


Figure  19 


The  right-handed  rectangular  coordinate  system  derives  its  name  from  the 
fact  that  a  right  threaded  screw  rotated  through  90®  from  OX  to  OY  will  advance  in  the 
positive  Z  direction,  as  illustrated  in  Figure  19.  The  system  is  right-handed  if  and  only 
if  the  rotation  needed  to  turn  the  x-axis  OX  into  the  y-axis  direction  OY  through  an  angle 
XOY  180  degrees  would  propel  a  right-handed  screw  toward  the  positive  side  of  the 
of  the  XY-plane  associated  with  the  positive  Z  axis  OZ. 

5.2.4.  Geocentric  Coordinate  System.  A  right-handed  system  with  its  origin  at 
the  center  of  the  reference  ellipsoid  of  the  earth,  the  positive  z  axis  in  the  direction  of 
geographic  north  pole,  the  x  and  y  axes  in  the  plane  of  the  equator,  the  positive  x  axis 
in  the  plane  of  the  meridian  through  Greenwich,  Figure  20. 

5.2.5  Local  Ground  Coordinate  System.  A  right-handed  system  with  the  x  and 
y  axes  in  a  plane  tangent  to  the  reference  ellipsoid  of  the  earth,  the  positive  y  axis  in 
the  meridian  plane  through  the  origin  and  pointed  to  the  North,  the  positive  z  axis  pointed 
toward  the  zenith  and  containing  the  exposure  station,  Figure  20. 


Local  Groand 


Positive 

Geoceatric 

X-Axis 


Figure  20 


5.3  Use  of  Matrix  Notation  in  Coordinate  Transformation. 

5.3.1  Definitions.  A  matrix  is  a  rectangular  array  of  numbers  or  functions  (later 
defined  as  elements)  of  the  form: 


where  m  represents  the  number  of  rows  in  the  array  and  n  the  number  of  columns.  This 
arrangement  is  defined  here  as  matrix  A.  If  m  =  n,  the  matrix  is  a  square  matrix.  A 
matrix  of  which  all  elements  except  the  diagonal  elements  are  equal  to  zero  is  called  a 
diagonal  matrix: 


A  unit  matrix  is  a  diagonal  matrix  whose  diagonal  elements  are  all  I's,  and  a  null  matrix 
is  one  in  which  all  of  the  elements  are  equal  to  zero. 


5.3.2  Addition  and  Subtraction.  The  sum  or  difference  of  two  matrices  is  a  new 
matrix  formed  by  adding  or  subtracting,  the  elements  which  lie  in  corresponding  positions, 
such  as 


or 
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By  using  the  letters  A,  B,  and  C  respectively  to  represent  the  above  matrices,  these 
operations  can  be  written  as 


A  +  B.C  or  A.B=:C 

5,3.3  Multiplication  of  Matrices,  The  product  of  two  matrices  is  defined  as  the 
new  matrix  formed  by  multiplying  the  corresponding  row  elements  of  the  first  matrix  by 
the  corresponding  column  elements  of  the  second  matrix  and  adding  the  results: 
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or  written  in  shorthand  form 
AxB  =  C 

In  order  that  matrix  multiplication  of  A  times  B  be  possible,  the  matrices  A  and  B  must  be 
conformable,  that  is,  the  number  of  rows  in  A  must  equal  the  number  of  columns  in  B.  If 
A  is  of  the  size  m  by  p,  that  is  m  rows  and  p  columns,  and  B  is  p  by  n,  then  C  is  of  size 
m  by  n.  To  multiply  a  matrix  A  by  a  number  k,  each  element  in  A  must  be  multiplied  by  k, 
or 


an 

^12 

a  13 

a^ik 

ai2k 

auk 

k 
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a  22 

3  23 

a2ik 

a  22  k 

323  k 

_a3i 

^32 

a  32 

! 

aaik 

a  32k 

a33k 
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5.3,4  Examples  of  Matrix  Notation  in  Linear  Transformation  Equations.  Matrix 
notations  are  commonly  used  to  represent  linear  transformation  equations,  involving 
translation  and  rotation  of  coordinate  axes.  The  following  is  a  set  of  such  equations  with 
matrix  notation  used  in  geodetic  coordinate  transformation.  Numbers  in  parenthesis  indi¬ 
cate  the  size  (mn)  required  for  the  matrices  to  be  conformable: 

(1)  X  =  x’  cos  a  “  y  sin  a 

y  =  x'  sin  a  +  y  cos  a 

z  =  z' 

Matrix  Notations: 


"^1 

cos  a  -  sin  a  0 

x’ 

y 

sin  a  cos  a  0 

y' 

z 

_0  0 

Z 

(3x1)  (3x3)  (3x1) 

by  letter  representation  for  each  matrix 
X  =  AX’ 

(2)  p  =  -I-  y'  cos  <J)  +  z’  sin  6 

0  =  Oo  -  y’  sin  <^  +  z’  cos  6 
R  =  +  x' 

Matrix  Notation: 


P 

■p7 

0 

cos  6 

sin  6 

x' 

0 

= 

Oo 

0  - 

sin  o 

cos  o 

y' 

_R_ 

i 

0 

0 

_zl_ 

(3x1)  (3x1)  (3x3)  (3x1) 


by  letter  representation  for  each  matrix 

P  =  P,  +  AX’ 

(3)  X’  =  (R-R,) 

y’  =  (P  -  P^)  cos  o  -  (Q  -  Qo)  sin6 
z'  =  (P  -  Pq)  sin  9  +  (0  -  Oq)  9 
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Matrix  Notation: 


x“ 

1^ 

0 

0 

1 

y' 

cos  <i>  -  sin  <j>  0 

O-Oo 

z  ’ 

sin  (f)  cos  <l>  0 

_R-R_^ 

(3x1)  (3x3)  (3x1) 

by  letter  representation  for  each  matrix 
x'=  AP  -  AP„ 

(4)  X  =  x'  +  Ak  (x  -  Xp)  +  A'k  (y  -  yp)  +  D'kc 

y  =  y'  +  Bk  (x  -  Xp)  +  B'k  (y  -  yp)  +  E'kc 
z  =  z'  +  ck  (x  -  Xp)  +  c'k  (y  -  yp)  +  F'kc 
Matrix  Notation: 


X 

x' 

A  A'  D~ 

X  -  Xp 

y 

= 

y' 

+k 

B  B'  E' 

y-Yp 

z 

— 

z' 

c  c 

_c 

(3x1)  (3x1)  (3x3)  (3x1) 

by  letter  representation  for  each  matrix 
X  =  x'  +  k 

(5)  Direction  Cosines  relative  to  XYZ  System 

X  =  AjX'  +  ^27' +  XjZ' 

Y  =  Mi'x'  +  Hjy'  +  H3Z' 

Z  = 

Matrix  Notation: 


"x' 

= 

Xj  Xj  Xj 

”x  ~ 

Y 

=: 

1^3 

y* 

_Z 

Jl  *’2  *’3_ 

z* 

(3x1)  (3x3)  (3x1)  or  by 

letter  representation  for  each  matrix 
X  =  \x\ 
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6.  PRACTICAL  APPLICATIONS 


6.1  Introductory  Statement.  Many  problems  of  analytical  triangulation  can  be  solved  by 
a  series  of  coordinate  transformations.  Vector  and  matrix  algebra  have  proved  to  be  most 
useful  for  this  purpose  when  the  work  is  to  be  done  by  electronic  computers.  Tdiis  section 
provides  a  useful  introduction  to  the  subject  of  coordinate  transformation. 

6.2  The  Aerial  Photograph.  An  aerial  photograph  is  always  considered  to  be  tilted. 

The  observed  rectangular  coordinates  of  the  image  on  a  tilted  photograph  are  called  x, 
y  and  z  with  z  as  the  focal  leng  .h  of  the  aerial  camera.  The  photograph  should  be  con¬ 
sidered  a  diapositive  with  the  x-  and  y-  axes  perpendicular  lines  defined  on  the  photo¬ 
graph  by  the  images  of  the  fiducial  marks  in  the  camera.  The  positive  x-  axis  is  the  axis 
which  most  closely  coincides  with  the  direction  of  the  flight  line.  The  y-  and  z-  axes  are 
chosen  to  make  a  right-handed  coordinate  system.  Thus  the  image-space  x,  y,  z  coordin- 
nates  of  the  nodal  point  are  0,  0,  -f  where  f  is  the  focal  length  of  the  lens. 

In  certain  instances  where  the  measurements  are  to  be  made  on  the  emulsion  side 
of  the  original  negative  (as  with  glass  plate  negatives)  the  positive  direction  of  the  y- 
axis  should  be  reversed.  The  image  space  x,  y,  z  coordinates  of  the  nodal  point  for  this 
case  are  0, 0,f. 

6.3  Three  Pairs  of  Conjugate  Axes.  In  addition  to  the  x,y,z  coordinate  system  which 
identifies  the  position  of  the  image  on  an  aerial  photograph,  an  auxiliary  system  (x’,y*,z’) 
is  introduced.  The  x*,y',z*  system  has  the  same  origin  as  the  x,y,z  system  and  differs 
from  it  only  in  orientation.  In  this  regard,  the  x',y',z'  axes  are  selected  parallel  to  the 
X,Y,Z  axes  of  the  local  (geodetic  or  geocentric)  ground  coordinate  system  where  +  Z  is 
upward  and  X  and  Y  form  a  right-handed  system  such  that  +X  is  East  and  +Y  is  North. 

The  angular  orientation  of  one  system  with  the  other  is  completely  specified  by  the  angles 
between  three  pairs  of  conjugate  axes. 


Figure  21 

Two  independent  systems  with  a  common  origin 
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Solid  Analytic  Geometry  of  the  Aerial  Photograph 


Figure  22 

Figure  22  shows  a  tilted  aerial  photograph  with  the  images  a,  b,  and  c  of  three 
ground  objects  A,  B,  and  C.  The  ground  objects  can  be  regarded  as  control  stations  for 
which  the  X-  and  Y-  coordinates  are  known  with  respect  to  some  oriented  geographic  sys¬ 
tem,  and  the  elevation  Z  is  also  known  relative  to  a  recognized  horizontal  plane  such  as 
sea  level.  The  geographic  coordinates  of  the  object  A  may  be  symbolized  as  X^,  Y^,  and 
Z^,  and  the  same  for  the  B  and  C.  The  perspective  center  0  of  the  photograph  is  also 
considered  to  have  geographic  coordinates  X^,  Yq,  Z^in  the  same  ground  system.  It  is 
convenient  in  most  problems  to  regard  0  as  the  origin  of  the  system  of  ground  coordinates. 

Figure  23  illustrates  the  geometric  relationship  between  a  spatial  object  A,  a 
photographic  image  B  as  recorded  on  an  image  plane,  and  a  camera  station  C.  In  the 
illustration  the  camera  station  (C)  is  the  center  of  projection  of  a  three  dimensional  space 
object  (A)  on  a  two  dimensional  image  plane  (B)  and  the  points  A,  B  and  C  are  considered 
to  be  on  a  straight  line. 
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6.4  The  Three  (Euler)  Rotation  Angles  co,  <j>  and  K.  The  angles  w,  <p  and  K  are 
defined  by  three  successiv^^  rotations  of  the  initial  set*  of  coordinates  x,y,z  into  an  aux¬ 
iliary  set  of  coordinates  x',y\z'.  The  positive  direction  of  the  rotations  is  considered 
counterclockwise  about  the  positive  direction  of  the  area. 

The  first  rotation,  omega  (<o),  or  x-  tilt  is  about  the  x-axis  (primary  axis)  con¬ 
sidered  counterclockwise  positive  as  viewed  from  the  plus  end  of  the  x-  axis  (Rosenfield 
considers  the  rotation  clockwise  as  viewed  from  the  origin  which  is  physically  identical). 
Omega  (to)  is  the  dihedral  angle  between  the  photo  perpendicular  and  the  XZ  plane.  That 
component  of  tilt  which  is  measured  around  the  x-  axis  is  called  roll.  Positive  roll  occurs 
when  the  right  wing  of  the  aircraft  is  lowered  in  flight. 

The  second  rotation  (Figure  24),  phi  (<|»),  or  y-  tilt  is  about  the  inclined  y-  axis 
(secondary  axis).  Phi  (6)  is  the  dihedral  angle  between  the  photo  perpendicular  and  the 
ZY  plane.  That  component  of  tilt  which  is  measured  around  the  y-  axes  is  called  pitch. 
Positive  pitch  occurs  when  the  nose  of  the  aircraft  is  lowered  during  flight. 

The  third  rotation.  Kappa  (K),  or  swing  is  the  rotation  about  the  twice  rotated  z* 
axis  (tertiary  axis)  which  is  the  perpendicular  lens  axis  or  the  z-  axis  of  the  photograph. 
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P  (Image  point) 

X 


Figure  24  The  <j)  -  rotation 

=  OM  =  00  -  MO 
=  00  -RM' 
z,,  =  MP  =  MR  +  RP 
=  0M*+RP 

OQ  =  OM’  cos  <|)  =  ^ 

QM’=  OM'  sin  <f>  - 
RP  =  M'P  cos  (j)  =  z^j^cos 
RM'=  M'P  sin  <f)  =Z^^sin4) 

y 4a3  ~  y^co</> 

Put  in  matrix  form,  the  last  three  equations  are: 


The  matrix  form  for  the  omega  and  kappa  rotations  are  derived  in  the  same  way.  These 
angles  can  most  conveniently  be  obtained  in  terms  of  the  three  direction  cosines  of  each 
of  the  three  axes  of  one  system  in  terms  of  the  other.  The  nine  cosines  are  sometimes 
spoken  of  as  the  vectorial  elements  of  space  orientation. 

6.5  Matrix  Form  of  Rotation  Formulae— Direction  Cosines.  The  rotation  formula  can 
readily  be  expressed  in  matrix  form.  (See  Figure  24.) 

The  matrix  for  the  primary  rotation  through  the  angle  omega  (<o)  about  the  x*axis 
is 


A(a))  = 


1  0  0 

0  cos  w  sin  w 

0  -  sin  CO  cos  CO 


This  rotation  generates  a  set  of  coordinates  which  is  then  transformed  into  another  set 
by  the  secondary  rotation  through  the  angle  phi  (^)  about  the  y-  axis. 

The  matrix  for  the  secondary  rotation  about  the  new  position  of  the  y-  axis  is 


AW 


cos  0  -sin  <^» 

0  1  0 

sin  0  cos  <t> 


This  rotation  generates  another  set  of  coordinates  which  is  then  transformed  into  still 
another  set  by  the  tertiary  rotation  through  the  angle  kappa  (K)  about  the  z-axis.  The 
matrix  for  the  tertiary  rotation  about  the  new  position  of  the  z-  axis  is 


A  (K) 


cos  K  sin  K  0 

-sin  K  cos  K  0 

0  0  1. 


The  three  successive  rotations  must,  of  course,  correspond  to  a  single  orthogonal 
transformation.  Since  the  sub-matrices  (A^o,  A<^,  AK)  refer  to  a  counter-clockwise  rotation 
of  a  coordinate  system  the  sequential  matrix  multiplication  to  form  a  single  matrix  must 
be  clockwise.  (AK,  A<^,  Aw)  as  follows: 


cos  K  sin  K  0 

-sin  K  cos  K  0 

0  0  1 


cos  <f>  0  -sin  <j> 

0  1  0 

sin  (f)  0  cos  <j> 


10  0 
0  cos  (0  sin  wl 

0  -sin  w  cos  CO 
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Applying  the  rule  of  matrix  multiplication  the  following  nine  elements  or  direction  cosines 
are  obtained: 


cos  <t>  cos  K  I  cos  w  sin  K 
-cos  <j)  sin  K  1  cos  co  cos  K 
sin  <f>  1 


+  sin  w  sin  <f>  cos  K  |  sin  w 

-  sin  w  sin  <f>  sin  K  |  sin  w 

-  sin  CO  cos  (f)  I 


sin  K  -  cos  CO  sin  <f)  cos  K 
cos  K  +  cos  CO  sin  <f}  sin  K 

cos  CO  cos  (f) 


6.6  Equations  Relating  the  Two  Systems  (x^  y,  z  and  x*,  y',  z*)  of  Rectangular  Coordi¬ 
nates  and  the  Nine  Direction  Cosines.  The  relations  between  the  nine  direction  cosines 
and  two  sets  of  rectangular  spatial  coordinates  (x,  y,  z  and  x',  y',  z’)  can  be  tabulated 
conveniently  in  the  following  conventional  manner: 


X' 

y' 

z ' 

X 

a  12 

^13 

y 

®  21 

3  22 

^  23 

z 

a  31 

3  32 

^33 

Exactly  the  same  information  is  given  by  conventional  vector  notation: 


^11 

^  12 

^13 

321 

^  22 

^23 

S  31 

^32 

^33 

The  inverse  notation  is  also  useful: 


X* 

^11 

321 

^31 

y 

CM 

03 

^22 

^32 

z'_ 

3  23 

^33_ 

(6.1) 


(6.2) 


Many  forms  are  used  in  mathematical  literature  to  represent  an  algebraic  matrix. 

In  this  text  a^n  represents  a  matrix  having  mn  elements,  arranged  in  m  rows  and  n  col¬ 
umns.  As  explained  in  section  5.3  the  first  subscript  indicates  the  row  in  which  each 
element  is  located  while  the  second  subscript  indicates  the  column  in  which  it  is  located. 
According  to  this  notation,  the  element  323  will  be  placed  in  the  second  row  of  the  third 
column  of  its  matrix.  The  range  of  the  subscripts  determining  the  number  of  rows  and  col¬ 
umns  defines  the  type  or  order  of  a  matrix. 
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The  elements  in  Equation  6.2  expressed  in  terms  of  the  direction  cosines  as  given  on  the 
top  of  page  40  are  as  follows: 

aji  =  cos  <l>  cos  K 

a  12  =  cos  (0  sin  K  +  sin  sin  ^  eos  K 

a  j3  =  sin  <0  sin  K  -  cos  uj  sin  ^  cos  K 

a2i  =  “  cos  ^  sin  K 

^22  =  cos  w  cos  K  -  sin  <o  sin  ^  sin  K 

a23  =  sin  (o  cos  K  +  cos  w  sin  <t>  sin  K 

a  31  =  sin  0 

332  =  -  sin  (0  cos  <l> 

£^33=  cos  03  cos 

The  Perspective  Transformation 


In  this  X'  y  z'-  coordinate  system  the  image  plane  is  considered  parallel  to  the 
XY-plane.  It  is  evident  from  the  similar  triangles  in  the  figure  that 

(X  -  X„)/(Z  -  Zo)  =  x'/z'  from  which  x'  =  (X  -  Xo)z'/(Z  -  Z,) 

Similarly,  by  visualizing  the  YZ-plane, 

(Y  -  Yo)/(Z  -  Zo)  =  y'/z'  from  which  y'  =  (Y  -  Yo)z'/(Z  -  Zo) 
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By  algebraic  identity 

(Z  -  Zo)/(Z  -  Zo)  =  z’/z’  from  which  z’  =  (Z  -  Zo)z’/(Z  -  Zo) 

By  substituting  these  values  for  x*,  y»  and  z'  in  Equation  6.2  and  dividing  the  expression 
for  X  and  y  by  that  for  z  we  obtain  the  following  basic  projective  transformation  equations: 

X  (X  ^  X  Ja,,  4-  (Y  ^  YJa,2  +  (Z  ^  Zja,3 

z  (X  -  X  Jaai  +  (Y  -  YJa32  +  (Z  -  Zjaaa 

(6.3) 

y  (X  ^  Xo)a2i  +  (Y  ^  YJa^z  +  (Z  ~  Z Ja23 

z  (X  -  X  Ja3i  +  (Y  -  Yo)a32  +  (Z  -  Zo)a33 

where  X,Y,Z  are  the  coordinates  of  an  object  on  the  ground,  Xq,  Yo,Zoare  the  coordinates 
of  the  camera  station  in  the  same  system  and  x,y,z  are  the  image-coordinates  in  which 
z  =  -  f  the  camera  focal  length.  Equation  6.3  expresses  the  condition  that  the  photographic 
image,  the  ground  object,  and  the  camera  lens  point  are  colinear.  Equation  6.3  is  transcen¬ 
dental  and,  in  the  most  general  case,  all  twelve  terms  are  considered  as  unknowns.  Con¬ 
sequently  a  form  of  Newton^s  method  is  used  to  solve  them.  This  is  an  iterative  method 
based  on  initial  approximations  which  are  quite  easily  obtained  for  all  the  unknowns.  The 
first  approximation  of  Phi  (^)  and  Omega  (<*))  is  always  zero.  Kappa  (K)  can  be  measured 
with  a  protractor  on  a  stapled  mosaic  or  a  photo-index.  The  flying  height  Z  is  known.  The 
X  and  Y  can  be  measured  with  sufficient  accuracy  for  the  problem  from  a  photo-index. 

After  these  substitutions  are  made,  the  computation  makes  all  the  other  corrections  by 
successive  approximations.  Usually  three  iterations  are  required. 

The  above  equations  are  too  involved  to  be  solved  for  each  of  the  six  parameters 
in  the  direct  manner  and  a  subroutine  computer  program  is  required. 

In  Equation  (6.1)  ai^,  ai2  and  denote  the  direction  cosines  between  the 
image-space  x-  axis  and  the  object-space  X  Y  Z  axis  respectively.  The  three  direction 
cosines  of  the  y-  axis  in  the  X  Y  Z  system  are  a 21 ,  Q22  and  323  and  similarly  for  the  z- 
axis.  The  three  direction  cosines  of  the  X-  axis  in  the  x  y  z  system  of  coordinates  are 
also  Uji  ,  a2i  and  a3i  and  similarly  for  the  other  axis. 

a^j  is  the  cosine  of  the  angle  between  the  X-  and  x-  axes,  322  is  the  cosine  of 
the  angle  between  the  Y-  and  y-  axes,  and  833  is  the  cosine  of  the  angle  between  the  Z- 
and  z-  axes.  Incidentally  333  is  the  cosine  of  the  tilt  angle  because  the  z-  axis  is  the 
photo  perpendicular  and  the  Z-  axis  is  the  plumb  line. 

The  nine  direction  cosines  of  the  matrix  are  not  totally  independent  values.  The 
sum  of  the  squares  of  the  elements  of  any  row  or  any  column  of  the  matrix  equals  one. 

The  sum  of  the  products  of  corresponding  elements  of  any  two  columns  or  any  two  rows 
is  zero,  or  more  specifically, 

^11  ^21  ^  ^12  ^22  ^13  ^23  =  0* 

The  value  of  the  determinant  of  the  matrix  is  either  +  1  or  -  1. 
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The  elements  of  the  matrix  are  the  elements  of  a  linear  orthogonal  transformation 
because  they  enable  the  computation  of  the  coordinates  of  a  point  in  either  system  from 
the  given  coordinates  of  the  same  point  in  the  other  system.  If  x,  y  and  z  are  known  for  a 
point,  and  if  the  nine  elements  of  the  matrix  are  known  then 

X  =  xaii  +  yaai  +  za^i 

or  if  X,  Y  and  Z  are  known  for  a  point,  and  if  the  nine  elements  of  the  matrix  are  known 
then 


X  =  Xaii  +  Yai2  +  Zai3 

6.7  Inverse  Case  —  Object  Space  Coordinate  System  in  Terms  of  the  Image  Space 
Coordinate  System.  The  matrix  expression  in  Equation  6.2  has  been  developed  for  the 
condition  wherein  the  image-space  coordinate  system  is  expressed  in  terms  of  the  object- 
space  coordinate  system.  The  coordinates  x,  y  and  z  can  always  be  measured  in  the  image 
coordinate  system.  If  the  inverse  case  is  desired,  the  object-space  coordinate  system  in 
terms  of  the  image-space  coordinate  system,  the  final  orientation  matrix  must  be  inverted. 
The  matrix  can  be  inverted  directly  by  transposition  of  its  elements  since,  for  an  orthog¬ 
onal  matrix,  the  inverse  is  equal  to  its  transpose.  The  inverse  form  can  also  be  developed 
by  inversion  of  each  of  the  individual  rotational  submatrices  which  would  necessitate 
inverting  the  order  of  matrix  multiplication. 

6.8  Relation  Between  the  Photographic  Coordinate  System  (Exposure  Station)  and 
the  Geocentric  Coordinate  System  (Local  Survey).  The  following  Equations  relate  the  two 
coordinate  systems: 


X 

an 

a  12 

a  13 

'x 

y 

_ 

a  21 

^  22 

a  23 

• 

Y 

z 

S  31 

S  32 

a  33 

Z 

where  x  y  z  are  photo  coordinates  X  Y  Z  are  geocentric  coor  Tnates 
X  =  longitude  (positive  in  value  east  from  Greenwich) 

<{>  =  latitude  (positive  in  value  north  of  the  equator) 


and  where  a. 


*22 


*23 


*31 


■*32 


=  sin  X 

COSX 

=  0 

=  -  sin</)  cos  X 
=  sin  f)  sin  X 
=  cos  ip 
=  cos  <j>  cos  X 
=  -  cos  <f>  sin  X 
=  sin  <l> 


(6.4) 
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In  the  geocentric  coordinate  system,  the  Z  axis  is  toward  the  north  pole,  the  +X 
and  -fY  axis  lie  in  the  plane  of  the  equator  with  the  +X  axis  in  the  meridian  plane  through 
Greenwich.  Therefore,  the  cosine  of  A  is  minus  for  longitudes  of  +90°  to  +180°  and  lati¬ 
tude  of  0°  to  +90°.  Consider  A  negative  in  the  western  hemisphere.  Formulas  for  this 
transformation  are  derived  by  rotating  and  translating  modified  geocentric  coordinates. 

The  classical  geocentric  coordinates  are: 

X  =  (N  +  h)  cos  sin  A 
Y  =  (N  +  h)  cos  ^  cos  A 
Z  =  [N(l  -  e^)  +  h]  sin<S) 

where  and  A  are  the  latitude  and  longitude  of  any  point,  h  is  the  elevation,  and 

N  is  the  length  of  the  normal  through  ,  A. 

Resection  in  photogrammetry  is  defined  as  the  determination  of  the  six  fundament- 
tal  parameters  w ,  (^) ,  K,  Xq  ,  Yq  ,  Z^  of  a  single  photograph  from  the  given  positions  and 
elevations  of  at  least  three  non-colinear  points  imaged  on  the  photograph.  The  X,  Y, 
Z“COordinates  in  Equation  (6.4)  are  in  the  form  required  for  resection.  The  final  adjusted 
X,  Y,  Z's  need  to  be  transformed  back  to  latitudes,  longitudes  and  elevations.  For  this 
the  following  formulae  are  needed: 

tan  A  =  X/Y 

tan  <(,  =  Z/(X^ 

Z  =  (N  +  h)  sin  ^ 

h  =  (X/cos  0  sin  A)  -  N  or 

h  =  (Y/cos  <f>  cos  A)  -  N  (use  the  equation  with  the 

larger  function  involving 
A) 

6.9  Orientation.  Since  orientation  is  most  important  in  this  analytic  system,  the 
following  brief  review  of  the  orientation  phase  is  included. 

The  orientation  phase  can  be  divided  into  interior  orientation  and  exterior  orien¬ 
tation.  Exterior  orientation  can  be  further  subdivided  into  relative  orientation  and  abso¬ 
lute  orientation. 

Interior  orientation  imposes  two  requirements: 

(1)  the  principal  distance  of  the  diapositive  must  equal  the  principal  distance  of  the  pro¬ 
jector,  which  is  the  length  of  the  perpendicular  from  the  interior  perspective  center  of  the 
projection  lens  to  the  plane  of  the  diapositive;  and 

(2)  this  perpendicular  must  intersect  the  plane  of  the  diapositive  at  the  principal  point  of 
the  photograph.  This  is  the  operation  known  as  centering  the  diapositive. 

Relative  orientation  is  defined  as  the  determination  of  the  three  angular  and  two 
linear  parameters  that  specify  the  attitude  and  position  of  the  photograph  (camera  station) 
with  respect  to  another  (overlapping)  one  that  shows  a  sufficiently  large  common  area. 
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Relative  orientation  is  perhaps  the  most  important  item  in  this  analytic  system:  it 
embodies  all  the  basic  mathematics  that  is  peculiar  to  the  system,  is  utilized  again  later 
in  resection  and  the  block  adjustment,  and  requires  the  second-largest  computer  effort. 

It  is  in  relative  orientation  that  the  principles  of  projective  geometry  are  applied.  A  clas¬ 
sic  geometric  rotation  of  the  axes  in  three  dimensions  is  needed  in  relative  orientation  to 
express  the  attitude  of  one  photograph  to  another.  Instead  of  using  the  three  angles 
between  the  respective  axes  as  in  analytic  geometry,  a  system  of  three  sequential  rota¬ 
tions  is  used,  the  primary  one  w  about  a  horizontal  x-axis,  the  secondary  one  <f}  about  the 
once  rotated  y-axis,  and  the  tertiary  one  K  about  the  camera  axis. 

Absolute  orientation  is  the  orientation  with  reference  to  a  vertical  and  a  horizon¬ 
tal  datum.  Absolute  orientation  procedure  comprises  three  operations:  (1)  bringing  the 
model  to  the  proper  size,  or  scaling  the  model;  (2)  leveling  or  horizontalizing  the  model; 
and  (3)  positioning  the  model  with  reference  to  the  horizontal  datum. 
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